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Abstract 

We present lattice simulation results corresponding to an SU{2) pure gauge theory 
defined on the orbifold space x Ji, where i?4 is the four-dimensional Euclidean space 
and /i is an interval, with the gauge symmetry broken to a f/(l) subgroup at the two ends 
of the interval by appropriate boundary conditions. We demonstrate that the U{1) gauge 
boson acquires a mass from a Higgs mechanism. The mechanism is driven by two of the 
extra-dimensional components of the five-dimensional gauge field which play respectively 
the role of the longitudinal component of the gauge boson and a massive real physical 
scalar, the Higgs particle. Despite the non-renormalizable nature of the theory, we observe 
only a mild cut-off dependence of the physical observables. We also show evidence that 
there is a region in the parameter space where the system behaves in a way consistent 
with dimensional reduction. 



1 Introduction 

Spontaneous symmetry breaking (SSB) is the phenomenon where the ground state of a 
system does not access all of its available symmetry, apparently breaking the symmetry 
group to a subgroup. In the Standard Model (SM) this is a crucial mechanism and it is 
not only responsible for predicting the existence of a fundamental scalar field, the Higgs 
particle, but also for the gauge bosons and fermions acquiring a mass. The somewhat 
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unsatisfactory fact about this mechanism in the SM, is that the Higgs potential, which 
is the concrete object that drives SSB, is input by hand at tree level in the Lagrangian, 
simply because we do not have any more fundamental way to generate it. There are 
many ideas of course trying to suggest an origin for the Higgs and its potential, one of the 
most elegant being that the Higgs field is the extra dimensional component of a higher 
dimensional gauge field and that the potential is generated quantum mechanically [1]. 
The earliest scenarios considered as extra-dimensional space the sphere S*^ [2-5]. In later 
applications the extra- dimensional space was taken to be non-simply connected, like or 
[6-8], so that the (non-contractible) Polyakov loops arc non-trivial. This is the general 
context where we would like to put ourselves in the present work. Alternative ways to 
achieve SSB in gauge theories with extra dimensions are discussed in [9,10]. 

There are enough motivations to take this idea seriously besides the economic way of 
generating the Higgs with its potential but the property that drew a lot of recent attention 
to these theories is the so claimed attractive possibility of all order finitcncss of the physical 
scalar mass. This sounds like a paradox from the beginning since the very point which 
has kept many field theorists rather hesitant from taking such an idea seriously is that 
higher (than four) dimensional gauge theories are non-renormalizablc and in a typical 
non-renormalizable theory one would expect that a mass parameter receives quantum 
corrections appearing in an arbitrary power of some dimensionless quantity built out of a 
dimensionful coupling and the cut-off. This is to be compared with the renormalizable SM 
where the couplings are dimensionless and the Higgs mass receives only a quadratic ultra- 
violet (UV) cut-off dependence under quantum corrections. Since even this quadratic 
UV sensitivity has been viewed as a drawback, supersymmetric generalizations of the SM 
were introduced and analyzed in detail, where the power like cut-off sensitivity is not 
present due to cancellations of infinities between superpartners. There is no doubt that 
supersymmetry is an elegant solution to the problem but it could happen that it is not 
realized at energies accessible in near future collider experiments so it is useful to be aware 
of alternative solutions. Back then to extra dimensions, in the case where the extra (fifth 
here) dimension is compactified on a circle, one can carry out a one-loop calculation of 
the Higgs mass and verify its aforementioned finiteness [11, 12] and can even give all order 
arguments to that effect [13-20], but the problem with this solution is that a simple circle 
compactification can not be realistic for various reasons, the absence of chiral fermions 
being one of the main. 

A way out is to compactify the extra dimension on an interval Ii instead of a circle, 
which can support chiral fermions at the two ends of the interval. Since the interval 
can be obtained easily from the circle by "orbifolding" , i.e. by identifying points and 
fields in the circle theory under the Z2 refiection operator TZ : — —x^, we will 
use the name orbifold when we refer to such a theory. A characteristic property of this 
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orbifold is that it is defined on a space with two four dimensional boundaries at each 
of the fixed points of the reflection action, where the gauge symmetry is reduced, thus 
naturally differentiating the boundaries from the rest of the space, which we call the bulk. 
An unfortunate consequence of field theories defined in such spaces is that the all order 
finiteness of the Higgs mass arguments are not anymore applicable because of the bulk- 
boundary interactions appearing at higher orders in perturbation theory which start to 
infect the finite bulk mass with cut-off dependence [21]. This is not unexpected; it is 
known that handling non-renormalizable theories analytically is not easy, in fact there is 
no general prescription that can be used in these theories such that their predictions are 
trustworthy. 

We would therefore like here to start a systematic investigation of higher dimensional 
orbifold gauge theories from the point of view of a lattice regularization [17,22,23]. The 
theory which will serve as our concrete example is an SU{2) gauge theory which has the 
symmetry broken (by boundary conditions) to its U{1) subgroup on the boundaries. This, 
we will argue, is a promising way of approaching extra dimensional theories: the goal is 
a non-perturbative understanding of a class of non-renormalizable theories and the hope 
is a non-perturbative understanding of the SM Higgs mechanism. 

It turns out that phenomenology puts surprisingly tight constraints on models and 
we will concentrate here on the two most immediate ones as one proceeds with the con- 
struction. One is associated with the mere existence of a four dimensional effective action 
and the second, closely related to SM phenomenology, the expected hierarchy of masses 
between the physical Higgs and the massive gauge bosons in the broken phase of the SM. 

1.1 Dimensional reduction 

The first issue is to show that an extra dimensional theory possesses in its parameter 
space a regime where it undergoes some kind of effective dimensional reduction. This can 
happen, in principle, in more than one ways. Historically, the first mechanism of hiding 
the extra dimension is going to the " Kaluza-Klein gauge" and making the size of the extra 
dimension small enough, so that its effects are negligible up to a certain energy scale and 
at the same time keeping the couplings perturbative so that the resulting effective theory 
can be useful for electroweak (or gravitational) physics. This is of course guaranteed as 
long as one treats the radius of the circle R, the cut-ofl0 A and the gauge coupling 
as independent parameters which is perhaps justified to a certain extent - it depends 
essentially on how far one is willing to go in perturbation theory. In a non-perturbative 
formulation on the other hand it becomes right away obvious that these three parameters 

^ Extra-dimensional theories are non-renormalizable and make sense only with a cut-off in place. 
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are tightly connected from the beginning and depending on where one sits in parameter 
space, a small change in one of the parameters could result in a dramatic change of the 
system. From this point of view, the region in parameter space where the theory describes 
the real world, (if it exists) could be only a small patch, which makes one wonder whether 
the free stretching of the dimensionful parameters sometimes employed to make a model 
phenomenologically viable is a valid operation. 

More recently, mechanisms of dimensional reduction which depend on localization 
rather than compactification have been proposed. In the context of gauge theories, all 
of these mechanisms involve one way or another a strong coupling and therefore non- 
perturbative physics. One, is the so called layered phase originally proposed in [24]. 
The idea is to investigate a five-dimensional lattice model where the gauge coupling in 
the four dimensional slices along the extra dimension is different from the gauge coupling 
that describes the interaction between the slices. This anisotropy could give rise to a new, 
so-called layered, phase where the static force is of Coulomb type in the four- dimensional 
slices and confining along the extra dimension, thus providing a localization mechanism. 
Some evidence for the existence of the layered phase in Abelian gauge theory was recently 
given in [25]. Another similar idea was due to [26] where it is assumed that the system 
has a phase where the bulk is in a confined while the boundary (defined by a domain wall) 
is in a deconfined phase, which forces the boundary gauge fields to remain localized. This 
idea was investigated on the lattice [27] and it was found that the low-energy effective 
theory contains not only the localized zero-modes but also higher Kaluza-Klein modes. 

A different mechanism of dimensional reduction was proposed in the context of the 
D-thcory rcgularization of non- Abelian gauge theories [28-30]. Here a non- Abelian gauge 
theory in five dimensions arises as a low-energy effective description of a five-dimensional 
quantum link model. The size L of four dimensions is taken to be infinity and the fifth 
dimension has a finite size R. If the five- dimensional theory happens to be in the Coulomb 
phase, the gluons are not massless but gain a mass which is exponentially small in the 
size of the extra-dimension. Thus by making R larger the correlation length given by 
the inverse gluon mass grows exponentially fast in R and therefore the extra-dimension 
disappears. 

In this paper we will not be able to give a conclusive answer to this important issue 
since it would involve a very extensive scan of the parameter space. We will give though 
some indirect evidence for dimensional reduction in the orbifold theory and we will leave 
the question of which mechanism is responsible for it, for a future work. At a qualitative 
level, we will be able to map the part of the parameter space of our model where sim- 
ulations were carried out onto the phase diagram of a known and well understood four 
dimensional theory, the Abelian Higgs model, which provides our first indirect evidence 
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for an effective dimensional reduction. Our more quantitative, but still indirect, evidence 
for dimensional reduction will be the measurement of the static potential between two 
infinitely heavy charged particles placed on four-dimensional slices located at the origin 
and in the middle of the fifth-dimension, separated by a distance r. Using the fact that we 
measure a massive U{1) gauge boson, we will fit the numerical data for the static potential 
to a five dimensional Yukawa potential and to a four dimensional Yukawa potential and 
compare the fits in both cases. 

1.2 Hierarchy of masses 

A universal feature that perturbative five dimensional orbifold pure gauge theories seem 
to posses is that the gauge bosons that survive on the boundaries are all massless and the 
mass of the Higgs turns out to be too low. In order to make some of the gauge bosons heavy 
one is forced to introduce fermions in the bulk [31], but still the ratio of the Higgs to the 
gauge boson mass is rrih/m^ « 1 [32,33] and thus phenomenologically excluded, unless 
additional assumptions are made. This property has been traced to the fact that since in 
five dimensions gauge invariance forbids a tree level Higgs mass, it has to be generated 
quantum mechanically. This one-loop mass turns out to be suppressed compared to the 
gauge boson mass which is governed by a vacuum expectation value of order one. More 
concretely, at one loop, the Higgs mass is found to be ruh ~ g^/R^^"^, where is the five- 
dimensional gauge coupling, whereas the gauge boson mass is basically a/R where 
a is a dimensionless vacuum expectation value obtained from minimizing the potential 
(see Appendix A). Realistic models require a very small a even though typical one- loop 
potentials tend to yield an a which is either zero or of order ~ 0.1. 

We will be able to compute the ratio uih/m^ non-perturbatively and find that, contrary 
to perturbative expectations, a Higgs mechanism is at work, resulting into a non-zero 
gauge boson mass, for a large part of the parameter space (in fact for the whole range 
of parameters we were able to scan). Again, disentangling the precise dependence of 
this ratio on compactification and finite lattice size effects requires a larger scan of the 
parameter space, which will be the topic of a future work. 

1.3 Organization of the paper 

In this paper we investigate numerically a one dimensional subspace of the parameter space 
of a five dimensional lattice gauge theory, starting from the vicinity of a first order phase 
transition and approaching the perturbative domain. In section [2], using generic properties 
of the measured orbifold spectrum, we argue that the spectrum of the system in part of 
this region seems to be consistent with an effective dimensional reduction. In section [3] we 
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construct in detail the lattice regulated theory and in particular its observables. In section 
IHwe present in detail our quantitative evidence for spontaneous symmetry breaking and 
dimensional reduction. In section [5] we state our conclusions. Finally, for completeness 
we provide two detailed appendices with some background, appendix El on 1-loop results 
and [B] on the derivation of the five- dimensional Yukawa potential. 

2 The orbifold effective theory 

We consider a pure SU{2) Yang-Mills (YM) theory in five dimensions with gauge potential 
Am, M = 0, 1,2,3,5. The fifth dimension is an interval, obtained by identifying points 
and fields on a circle of radius R under the Z2 reflection operator TZ : — > —x^. The 
projection breaks the gauge symmetry at the two ends of the interval according to 



In the picture of dimensional reduction as it is known in finite temperature field theory [34] 
the five-dimensional fields are expanded in a Fourier series in the quantized momentum 
along the compact dimension. At low energies, in our case much below the compactifica- 
tion scale given by the inverse radius 1/R, we have an effective four- dimensional theory 
of the zero-modes (i.e. constant along the fifth dimension) of the fields. 

But the Fourier series that defines the Kaluza-Klein expansion breaks gauge invari- 
ance. This we cannot afford at the non-perturbative level. Here the particle spectrum is 
read off from correlations of five- dimensional, gauge invariant operators. These operators 
are classified according to specific symmetries. If dimensional reduction occurs, a mass gap 
in the various particle channels should be seen. The ground state of the operators which 
corresponds to the gauge field component^] has the quantum numbers of a complex 
scalar and the ground state of the A^, = 0, 1, 2, 3 component^ have quantum numbers of 
a f/(l) gauge field from the point of view of four dimensions. Clearly, the lowest lying spec- 
trum of the orbifold theory should coincide with the one of the four-dimensional Abelian 
Higgs models, if dimensional reduction works like in finite temperature field theory. The 
spectra with the higher excitations included will be different, as in the five- dimensional 
theory they are sensitive to the compactification scale 1/R. Before defining the orbifold 
theory in great detail and present extensive simulation results (which will be the topic of 
the next sections), one would hope to be able to exploit this similarity of spectra by map- 
ping the phase diagram of the orbifold simulation on the phase diagram of the Abelian 

^ In Section [3. 1.21 we will specify the gauge-invariant meaning of this, related to the Polyakov line. 
In Section 13.21 we will construct a gauge-invariant operator that corresponds to the gauge bosons, 
very much in analogy with the four-dimensional SU{2) Higgs model [35]. 
^ We thank U.-J. Wiese for his comments in this respect. 
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Higgs model. Being able to do so, would be a first circumstantial evidence for an effective 
dimensional reduction. 

Next, we define the parameter space of our model. In a lattice regularization, the cut- 
off is provided by the inverse lattice spacing A = and it preserves gauge invariance. 
That is why the lattice formulation is particularly useful to study non-renormalizable 
theories, which requires a cut-off. The parameter space of the lattice theory is two di- 
mensiona|§ . One of the parameters is the (dimensionless) number of lattice points in the 
fifth dimension 

iV5 = ^ (2.2) 
a 

and the second is the dimensionless five dimensional lattice coupling 

^ 2N 2N , , 

(3 = -^a=-^, (2.3) 

i/5 i/O 

where go is the dimensionless bare gauge coupling. For later reference, notice that in 
regions of the parameter space where has a mild cut-off dependence, j3 is inversely 
proportional to A. All observables are computed in units of the lattice spacing, or equiv- 
alently in units of R, and are thus functions of and (3. One can also define a derived 
parameter, the coupling 

- 3 = f . (2-4) 

for the effective four-dimensional theory. It is also possible in principle to use an anisotropic 
lattice [36] in which case a lattice gauge coupling (and a different lattice spacing) is de- 
fined for the fifth dimension {(3^) and a different one for the four dimensional subspace 
(Pi). Here we will restrict ourselves to isotropic lattices for which ^4 = ^5 = (3. 

From simulations of the five dimensional SU{2) gauge theory on the orbifold 3^/7^2 
with the number of points in the extra dimension fixed to = 4, and on an isotropic 
lattice with an associated coupling /?, the spectrum can be safely determined when j3 
is larger than a critical value jSc = 1.5975, which separates a confined {(3 < jSc) from a 
deconfined (/5 > Pc) phase. The spectrum measured in simulations corresponding to the 
deconfined phase consists of a massive Higgs ground state and a massive Abelian gauge 
boson, along with their excitations [23] which implies that our system in this region of 
parameter space crosses from a confined phase into a Higgs phase. □ In order to be able to 
make a comparison with the Abelian Higgs model we now recall a few well known facts. 



^ For this discussion we assume the four-dimensional sizes of the lattice to be infinite. 

^If the four-dimensional effective theory were just a pure U{1) gauge theory, (due to the orbifold 
breaking of the SU(2) symmetry) and dimensional reduction occurs, one should be able to map our 
results onto the four-dimensional Abelian model. For A's = 4 we have = /3. The pure compact U{1) 
gauge theory has a phase transition at f3c — 1-01 [37] and therefore we would end up in its Coulomb 
phase. Clearly, this can not be sufficient to describe our system because, as already mentioned, in the 
deconfined phase we measure not only massive scalars but also a massive gauge boson. 
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Figure 1: Phase diagram of the Abehan Higgs model for the Higgs field with charge q = 2. 
From studies at A = oo [38]. 

We summarize results from analytic [38, 39] and numerical [40] studies of the four 
dimensional Abelian Higgs model. The Euclidean action may be written as 

S = ^|a(01'(x)0(x) - 1)2 + 0^ (x)0(x) 

-K J2 {(pHx)e"'^-^''U{x + ajl) + c.c.) \-^J2^p- (^.S) 

A* P 

Here, is the complex Higgs field, q is its charge and k, the hopping parameter related 
to the bare mass mo through 

1 - 2A „ 



2 2 

a rriQ 



K 



(2.6) 



The analytic study in [38] was performed in the limit A ^ oo where the length of the 
Higgs field p is frozen to 1. In the unitary gauge the action then reads 



S = -«:5^2cos[gA,(x)]-|^f/,. (2.7) 

There are two interesting cases. One is when the Higgs field is in the fundamental rep- 
resentation, that is it has charge q = 1 and the second is when it has charge q = 2. For 
us it is the second case that is relevant since, as will show in the next section, the lattice 
operator which is identified with the Higgs has charge two. 

We inspect the action for A = oo in the unitary gauge, Eq. (12.71) . in the limit k ^ oo. 
In this limit the gauge variable A^{x) can take, for g = 2, two values, or vr, corresponding 
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to Z2 gauge links U{x,fi) = ±1. (For general charge q it will be a Zg gauge theory). The 
Z2 gauge theory has a second order phase transition. The analytic study of [38] in the 
A = 00 case concludes that there are three distinct phases, sketched in Fig. [H The main 
difference compared to the q = 1 case is a phase boundary that separates the Higgs from 
the confinement phase. The static potential in the Higgs phase is of Yukawa type. In the 
confinement phase it rises linearly (area law for Wilson loops) which is not the case in the 
confinement region for a fundamental Higgs. 

Clearly, this is the version of the Abelian Higgs model on which the orbifold theory 
can be naturally mapped since the lowest lying spectra are the same and the part of the 
orbifold phase diagram we have investigated can be recognized inside the phase diagram 
of the q = 2 Abelian Higgs model. 

3 Orbifold on the lattice 

Gauge theories on the orbifold can be discretized on the lattice [17,22]. One starts with 
a gauge theory formulated on a five- dimensional torus with lattice spacing a and periodic 
boundary conditions in all directions M = 0, 1, 2, 3, 5. The spatial directions (M = 1, 2, 3) 
have length L, the time-like direction (M = 0) has length T, and the extra dimension 
(M = 5) has length 27iR. The coordinates of the points are labelled by integers n = {um} 
and the gauge field is the set of link variables {U{n, M) G SU{N)}. The latter are related 
to a gauge potential Am in the Lie algebra of SU{N) by U{n,M) = exp{aAM{n)}. 
Embedding the orbifold action in the gauge field on the lattice amounts to imposing on 
the links the Z2 projection 

(l-r)f/(n,M) = 0, (3.1) 

where F = TZTg. Here, TZ is the reflection operator that acts asTZn = {n^, —n^) = n [fi = 
0, 1, 2, 3) on the lattice and as TZ U{n, jj,) = U{n, jj,) and TZ U [n, 5) = U\n — 5, 5) on the 
links. The group conjugation Tg acts only on the links, as TgU{n,M) = gU{n, M)g~^ , 
where g is a constant SU{N) matrix with the property that g^ is an element of the 
centre of SU{N). For SU(2) we will take g = —ia^. Only gauge transformations {f2(n)} 
satisfying (1 — T)VL = are consistent with Eq. (13.11) . This means that at the orbifold 
fixed points, for which = or = N^, the gauge group is broken to the subgroup 
that commutes with g. For SU{2) this is the U{1) subgroup parametrized by exp(z</)cr^), 
where are compact phases. 

After the projection in Eq. (13. ip . the fundamental domain is the strip Ii = {nf^,0 < 
"^5 < The gauge- field action on Ji is taken to be the Wilson action 

p 
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where the sum runs over all oriented plaquettes U{p) in /i. The weight w{p) is 1/2 ii p is 
a plaquette in the {^u) planes at ns = and = N^, and 1 in all other cases. Dirichlet 
boundary conditions are imposed on the gauge links 

U{n,n) = gU{n,iJ,)g~^ at = and = N^. (3.3) 

The gauge variables at the boundaries are not fixed but are restricted to the subgroup of 
SU{N), invariant under Tg. The Wilson action together with these boundary conditions 
reproduce the correct naive continuum gauge action and boundary conditions on the 
components of the five-dimensional gauge potential [17]. For example, for SU{2), 
("photon") and Ag'^ ("Higgs") satisfy Neumann boundary conditions and Aj^'^ and 
Dirichlet ones. 

One of the main results of [17] was to show, through a geometrical construction, that 
the orbifold projection equation Eq. (13. ip implies the absence of a boundary counterterm 
for the Higgs mass. Given the explicit breaking of the gauge invariance at the boundaries, a 
boundary mass term for A^ is invariant under the unbroken gauge group. In the continuum 
this mass term would be tr{[A5, (7] [A5, (7^^]} evaluated at the boundaries X5 = and 
X5 = ttR. If present, such a term would imply a quadratic sensitivity of the Higgs mass to 
the cut-off. For the lattice action Eq. (13.21) this would require to add a boundary action 
term with an additional coefficient ft. As the lattice spacing a changes, a fine tuning of fl 
would be required to keep the Higgs mass finite. Fortunately this term is absent and the 
orbifold action is simply Eq. (13.21) . Since five- dimensional theories are non-renormalizable, 
they make sense only as effective theories for energy scales much below the cut-off. On 
the lattice this is the Symanzik effective action [41-44], a continuum action which is 
a systematic expansion in the lattice spacing a. For the orbifold theory the Symanzik 
effective action is [17] 
1 



S 



2gl 



J d'ztl{FMNFMN} 

+abi d'^ X Re ti {g Fm N Fm n} + ab2 d^xRetT{gFMNgFMN} 

Jz={0,'rrR} J z=iO.TTR\ 

Wc J d^zil{DLFMNDLFMN] + • • • 



2={0,7r_R} J z={Q,-kR} 

(3.4) 



where Fmn is the field strength tensor, its covariant derivative and the coefficients 
61, 62, c, . . . are computable in perturbation theory. At 1-loop, &i 7^ and 62 = [11]. A 
boundary mass term for the Higgs would appear with a coefficient jl/a?. 

3.1 Operators for the Higgs 

If the fifth dimension were infinite, the gauge links U{n., 5) would be gauge-equivalent 
to the identity, which corresponds to the continuum axial gauge ^5 = 0. On the circle 
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Figure 2: The Polyakov line P on /'L2- 



S*^ one can gauge-transform U{n, 5) to an ns-independent matrix V{n,) that satisfies 
P = where P = P{n^) is the Polyakov line winding around the extra dimension 

at four- dimensional location n^. Therefore an extra-dimensional potential {A^)\a.t can be 
defined on the lattice, through V = exp{a(y45)iat}, as 

a{A,U = ^^{P-P^) + 0{a'). (3.5) 

At finite lattice spacing the O(a^) corrections in Eq. fl3.5p are neglected. One easily 
checks that TIP = P"^ and so 7^ (y45)iat = — (^5)iat5 as it should be to have the same 
transformation behaviour as in the continuum. 

In order to construct the gauge potential A^ on the orbifold 5'^/Z2 we start from the 
circle parametrized by the coordinates = —N^, . . . , A^s — 1 [N^ = irR/a is identified 
with — A^s). We impose on the links building the Polyakov line P 

Pin,) = Uiin,, 0), 5) . . . ?7((n^, N^-l), 5) ?7((n^, -A^5), 5) . . . -1), 5) (3.6) 

the orbifold projection Eq. (13.11) . The result is 

P{n,) = Kn,)gl\n,)g-\ (3.7) 

with 

l{n,) = ?7((n^,0),5)f/((n^,l),5) ... f/((n^,Ar5-l),5). (3.8) 

The Polyakov line Eq. (13. 7p on /'L2 is shown schematically in Fig. [2] and from it we 
define {A^)\sx{n,) = {P{n,) — P'l'(n^)}/(4A'5). Being anti-Hermitian, the field {A^)i!^t can 
be represented using the unit matrix and the Hermitian generators of SU{N): 

(/l5)iat = -tgo{AllN + AfT^). (3.9) 

In order to construct the Higgs field on S^/Z2 we have to project^ (A5)iat onto the com- 
ponents A^ for which gT"- g^^ = — (the other generators have gT"" g~^ = T"). The 

^ In the continuum the orbifold projection selects automatically the components . On the lattice, 
due to the lattice artifacts 0{a^) in the definition Eq. p.Sp . "wrong" components {^5)1^,1 can be non-zero. 
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Higgs field is now defined as in the continuum 

$(n^) = [a{A,),,i{n^),g] = 2iagoAl{n^)gT\ (3.10) 

Note that the commutator projects out the identity component of (^5)1^1- Under gauge 
transformation Q{n), l{n^) — > f2(n^, 0) /(n^) f2(n^, A^s)""^. Since [^l{n),g] = at the 
orbifold boundaries rig = and — N^, it follows that 

n(n^,0)$K)(](n^,0)-^ (3.11) 

The Higgs field $ transforms like a field strength tensor at the boundary. In the special 
case of gauge group SU(2) only a [/(I) gauge symmetry survives at the boundaries. If 
we parameterize the U{1) boundary gauge transformations by 0(n^, 0) = exp{ia;(n^)cr^}, 
the Higgs field transforms as 

«=f , ° = \ - f " ^"M. (3.12) 

showing that it has charge 2 under the U{1) gauge group. 

Since tr{$} = 0, in order to extract the Higgs mass we define 

H{no) = (f)' E tr{$(n,)$t(n^)} (3.13) 

ni,n2,n3 

and build the connected correlation 

C{t) = {(^(^0) ^(^0 + ^/«)) - (^(^0)) (^(^0 + ^/«))> (3.14) 

no 

const. X e^'"''* , 

and the Higgs mass can be extracted from the effective masses amh^es{t + a/2) — 
\n{C(t)/C(t + a)}. Writing the correlation C(t) as 

C{t) = {{H{t)-{H)){H{0)-{H))) (3.15) 

one can see that it is a sum of positive and negative numbers of order one, the result 
being a small number due to cancellations. On the other hand, the variation 

AC = {{H{t)-{H)f{H{0)-{H)f)-C^ (3.16) 

is a sum of positive numbers of order one minus a very small number and hence of order 
one. Furthermore, AC is essentially independent of t. This means that the error in the 
effective Higgs masses is approximately 



Aa™.,e. . AC J(^) + (^) , (3.17) 
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and, since C(t) ~ e its t-dependence is 

Aamh,es ~ e"^'^*. (3.18) 

In the above, rrih is the plateau value of the Higgs mass which is constant and therefore 
one expects that the error in rrih increases with t exponentially. 

3.1.1 Variational technique 

It turns out that the correlation function Eq. (13.141) suffers from a loss of significance in 
numerical simulation. This is the, unfortunately, common problem of signals, which are 
exponentially small in the time t, with an almost constant variance, and hence constant 
statistical error. The ratio of the signal to the error falls off exponentially in the time t; in 
the large-t region, where the leading exponential decay Eq. (13.151) due to the Higgs mass 
should be seen, the signal is lost. This is reflected in the exponentially growing error of 
the effective masses Eq. (I3.18p . 

To cure this problem we employ the variational technique of [45]. A basis of Higgs 
operators is constructed and the best operator to capture the Higgs mass in Eq. (13.151) . 
i.e. whose overlap with the eigenstate of the Hamiltonian is the largest, will be a linear 
combination of these basis fields. If the contributions from excited states are suppressed 
then the leading exponential Eq. (13.151) can be extracted at smaller values of t, where the 
signal might not be lost. 

Here we sketch how this works; more details can be found in [45]. We construct a 
set of Euclidean fields Oi, i = 1, . . . , r with the same quantum numbers as H = Oi in 
Eq. (I3.13p . Then we build the matrix correlation function 

Q,{t) = {o.{t)o,m - mt)) {o,m , (3.19) 

and write the spectral decompositions 

mt)0,m = ^^e-^"^-*^^'"-^")^^^.* (3.20) 

m,n 

mt)) = ^5^e-^"^A», (3.21) 

n 

where 

Z = 5^e-^"^ and A» = (m|O,(0)|n). (3.22) 

n 

Here m,n = 0, 1, 2, ... label the eigenstates \m) with energy eigenvalue Em of the Hamil- 
tonian EI and T is the temporal size of the lattice. We use the same symbol Oi to denote 
the Euclidean field and the corresponding operator in the Hamiltonian formulation. 
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There are two effects, which derive from the finiteness of T and the periodic boundary 
conditions in time (which imply taking the overall trace in the spectral decomposition) 
[46]. Firstly, if the operators Oi have a non- vanishing expectation value the connected 
correlation functions Eq. f l3.19p have in general t-independent contributions 



^00 ^11 



aH) _ aU) 
^00 ^11 



(3.23) 



where nih = Ei — Eq is the mass gap, i.e. the Higgs ground state mass. 

Secondly, in the limit that T, t are both large and the difference T — t is close to t, 
the matrix correlation function Eq. (13.191) has the leading behavior 

C^,{t) ^ 5^ [e-(^"-^°)* + e-(^"-^")(^-*)] ^2^;^^ + 

n>0 

^ ^-iEr.-Eo)T-[E^-E^)t^(r^^^i,)* (3^24) 
•m,n>0 

Here we have assumed that Amn = A*^m ^'^^ Cij(t) is real, which is true if the operators 
Oi are Hermitian. If we take t = T/2, the contribution of the second term in Eq. (13.241) 
goes like exp{ — [{En + Em)/2 — Eo]T} compared to exp{ — (_E„ — _Eo)T/2} of the first term. 
The second term is hence subleading and can be neglected. We get 



C^At) — Cl^it) + Cl^iT - t) , C:,(t) = $^A2A2V(^"-^«)*. (3.25) 

n>0 



The correlation function is a sum of two contributions and is symmetric about t = T/2. 

First we assume that T is large enough so that the t-independent contribution Eq. (I3.23P 
is negligible and that t is small enough so that the contribution T — t is also negligible. 
Then the correlations Eq. (13.191) have the T = oo behavioi^ (assumed in [45]) 

oo 

a,{t) = ^A«aSV*^% W. = E^-Eo. (3.26) 

a=l 

The lowest masses Wa can be extracted by solving the generalized eigenvalue problem 

Cit)i,^lJaAt,to) = A„(t,to)ai(to)V^a,,(t,to), (3.27) 

where the correlation matrix is taken at variable time t and at fixed time to (we may set 
to = 0). From a mathematical lemma proved in [45] it follows that 

A,(t, to) *=°°c,e-*^<^ [1 + 0(6-*^'^'^)] , a = l,...,r, (3.28) 



Then in Eq. ([X^ and Eq. ([X^ only the term n = contributes. 
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where > and AWa = min^^Q, \ Wa — W/jj. One expects that Cq ~ QtoWa g^^^ 
the coefficients of the correction terms in Eq. (I3.28p . because of the excited states, are 
suppressed. The masses can be extracted from 

at moderately large value of t. 

If the contribution T — t is not negligible, as it happens when t approaches T/2 (and 
we might be forced to go to such large values of t to find a plateau), then a different 
formula has to be used. The starting point is Eq. (13.251) . If we say that A^(t, to) solves the 
generalized eigenvalue problem Eq. (13.27^ for the matrix C defined in Eq. (13.25^ . then it 
is straightforward to show that 

Xa{t,to) = X'^{t,to) + X'^{T -t,to) (3.30) 

solves the generalized eigenvalue problem for the full matrix C. Using Eq. (I3.28P for A' 
we get the formula 

Xa{t,to) = 2cae-^-^/\osh[{T/2-t)Wa]. (3.31) 
By using the ratio ri2 = Xaiti^to) / Xaih^to) and the definitions 

X = e-^- , ri = T/2 - ti , rs = T/2 - , (3.32) 

it is easy to arrive at the equation 

ri2 (x^^ + x-^') - {x^' + x-^') = (3.33) 

to be solved numerically (Newton-Raphson) for x. We take = 0, ti = t and t2 = t + a 
for t = a, 2a, . . . , T/2 — a. 



3.1.2 Higgs operators 

A basis of Higgs operators Oi = Hi, defined as in Eq. (I3.13p . can be constructed by 
modifying the definition of the Higgs field $ to create a set of fields We can for 
example consider displaced Polyakov lines on S'^/Za as shown in Fig. [3l The position 
where the displacement in one of the spatial directions k = 1,2,3 takes place can be 
varied along the extra coordinate = 0,1, . . . , \_N^/2\. The displacements at the other 
ns values are equivalent by symmetry with respect to reflection — > — n^. For the 
displacement at ns = Eq. (13. 7p is replaced by 

P^'^M = I J2 Uin,,n'^)\^^=oWU^in,,n'^)U,=^.^gl\n,)g-\ (3.34) 

\n—n'\=a 
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and similar expressions P("5) forns = 1, . . . , [N,J2\ . In Eq. flS^Ij) U{n^, n'^) is the parallel 
transporter from the four-dimensional point n'^ to in the indicated slice in the extra 
dimension. 

Yet another possibility to create fields for the variational basis is to consider smeared 
Higgs fields. The simplest is shown in Fig. HI the field (represented by a thick 

point) is replaced by a smeared field $'(72^) made of a linear combination of and 
the nearest neighbour fields in the three dimensional space 



^\n,) = (l-a)^n,) + - f^K.<)k=o f/^K,Olns=o • (3.35) 



The smearing parameter is a. The definition of Eq. (13.351) ensures that $' transforms as 
$ under the gauge transformation in Eq. (13.111) . Variants of the smearing technique can 
be found in [47] . The smearing procedure Eq. (I3.35p can be iterated a number of times, 
which we take to be 3. We set a = 0.7. 

Finally, the gauge links used to construct all the Higgs fields are replaced by smeared 
gauge links. This is a very simple procedure to create extended operators. We have 
implemented APE smearing [48] for the spatial links, i.e. U{z, k), k = 1,2,3 and U{z, 5). 
The links are decorated with staples in the spatial directions / = 1,2,3, clearly with the 
restriction / 7^ k for U{z,k). The smeared link U' is obtained by adding to (1 — a)U 
the sum of the decorating staples multiplied by a/ (number of staples), where a is the 
smearing parameter. The smeared link is projected back onto SU{2). The smearing of 
the gauge links is iterated 3 times with a = 0.75. 



"0="o 
n—n'\=a 





1 




+ 







Figure 3: The displaced Polyakov line L on 3^/7^2. 
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Figure 4: A smearing procedure for the Higgs field $. 
3.2 Photon operators 

In this section we describe the construction of operators in the SU{2) orbifold, which 
create the gauge boson (photon) associated with the unbroken U{1) gauge group on the 
boundaries. For each Higgs field $j constructed as described in Section 13.1.21 we define 
the SU (2) valued quantity (in the following we suppress the index i) 

«W = //wl'? (3.36) 
v/det($(n^)) 

We note that from Eq. (13.101) and g"^ = —g (which holds for SU(2)) it follows = —a. 
For the three spatial directions = 1, 2, 3 we define a "decorated double link" 

V{n^,k) = U{n^, k) a{n^ + ak) U\n^, k) a{n^) (no sum over fc) , (3.37) 

which transforms like $ under gauge transformations. In analogy with the definition [35] 
in the standard SU (2) Higgs model, we define an operator for the photon field 

Wkirif,) = -itT{a^V{nf,,k)}. (3.38) 

We will consider this field projected to zero three-dimensional momentum 

Wk{no) = (I)' Yl ^fcK)- (3-39) 

It is not difficult to show that the operator H4(^o) is real, invariant under the group 
conjugation Tg and odd under the three-dimensional parity transformation. In order to 
study the naive continuum limit, we define a continuum gauge potential at the boundaries 
through 

U{x,k) = e'^^'=(^) = l2 + aAfe(x) + 0(a2), (3.40) 

where Ak = —igoAla^. Here we use the continuum notation to label the lattice points, 
Xfj^ = n^a. Using the definition of lattice derivative a{x + ak) = a{x) + adka{x) and the 
properties a{x)^ = — 12 and {Ak{x),a{x)} = we get up to O(a^) corrections 

Wk{x) = zatiiasaix) [dk + 2Ak{x)] a{x)} . (3.41) 

Dk 
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Consistently with Eq. (13.121) the covariant derivative Dk has a charge 2 in front of Ak. In 
summary, in the naive continuum hmit Wk{x) corresponds to a covariant derivative term 
for the Higgs field and is hence a vector in the representation of spin 1. 

The photon mass can be extracted as follows: for each Higgs field in the varia- 
tional basis discussed above we construct using Eq. (13.381) the corresponding photon field 
W''^\nQ). Then we build the matrix of connected correlation functions 



The expectation value of {no) is zero on average, nevertheless it is subtracted to 
possibly reduce the statistical noise. The matrix {t) is treated as its counterpart for 
the Higgs to finally get the masses. 

3.3 The static quark potential 

We measure the potential between a static quark and a static antiquark placed in the 
four-dimensional slices as we vary the location = 0, 1, ... , \_N^/2\ of the slice. The 
potential is extracted from the expectation values of the traces of rectangular Wilson 
loops W{t,r) of size r in three-dimensional space and t in time direction. The Wilson 
loops are averaged over the three spatial directions and smeared gauge links, constructed 
like it is explained in Section 13.1.21 are used for the spatial Wilson lines. The potential 
V{r) is then defined as the plateau value at large time t of the effective masses 



In the boundary slices of the orbifold the gauge links belong to gauge group U{1), whereas 
in the bulk slices they are SU{2) gauge links. Since the links on the boundary are "colder" 
(i.e. the boundary four- dimensional plaquettes have a larger value) [22], the potential 
turns out to be more precise on the boundary. Also, as the value of is increased, the 
statistical errors on the potential also increase. This calls for improving the extraction of 
the potential and there are methods to do this. 

4 Numerical results 

We present detailed simulation results of the SU{2) gauge theory on the orbifold. The 
main part of the simulations is a scan in /3 at fixed = 4. The three-dimensional 




(3.42) 




(3.43) 
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Figure 5: The phase transition on the orbifold at = 4. The plot on the left hand side 
shows the behaviour of vacuum expectations values (tr{$''"$}) for different Higgs fields. 
The plot on the right hand side shows a met ast ability in the Polyakov line in one of the 
three-dimensional directions right at the phase transition /3c = 1.5975. 



space size is L/a = 8 and the temporal size is T/a = 96. The large value of the latter 
was necessary in order to extract reliably the Higgs mass, which turns out to be very 
close to the 1-loop perturbative value. To check for finite L effects we also performed 
for some (3 values simulations at L/a = 12. The algorithm is based on SU{2) heatbath 
and overrelaxation updates in the bulk and U(l) heatbath and overrelaxation updates on 
the boundaries. The statistics of the simulations at L/a = 8 varies between 90'000 and 
260'000 measurements of the observables, at L/a = 12 it is of 32'000 measurements. Each 
measurement is separated by one heatbath and L/2a overrelaxation sweeps (i.e. updates 
of all the gauge links). 

We have also run some simulations to measure the static potential in the four-dimensional 
slices along the extra dimension. For these runs the orbifold geometry was chosen to be 
T/a = 32 and L/a = 16 and we compare the potentials with = 4 and A^s = 6 at 
f3 = 1.609. This f3 value was chosen, since the mass spectrum there resembles more what 
we expect in a compactification scenario. 



4.1 Phase transition 



In infinite volume or with periodic boundary conditions on a finite torus five- dimensional 
SU{N) gauge theories have at least two phases [36,49,50], a confinement massive phase 
at small values of (3 and a deconfinement or Coulomb massless phase at large values of jS. 
For SU{2) the phase transition is located at Pc = 1.64 [36,49]. With orbifold boundary 
conditions this phase transition persists but the critical value jSc depends on and on 
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L/a. It is signalled by a jump in the expectation values of plaquettes, see [22]. In Fig. [5] 
the plot on the left hand side shows what happens on the orbifold with A^5 = 4 and 
L/a = 8 with vacuum expectation values tr{$$^}, where for $ we take some of the Higgs 
fields as explained in the legend and in Section [3l There is a clear discontinuitjQ in the 
vacuum expectation values at = 1.5975. 

There is a strong indication that the phase transition is of first order. The plot on 
the right hand side of Fig. [5] shows the historj0 of the expectation value of the Polyakov 
line in one of the three-dimensional space directions evaluated at the boundary = of 
the orbifold. As it was shown in early Monte Carlo study of SU{2) gauge theory at finite 
temperature [51,52] the expectation value of the Polyakov line is zero in the confined phase 
and it becomes nonzero in the deconfinement phase, thereby breaking spontaneously a 
Z2 symmetry which changed the sign of the Polyakov line. In finite volume one actually 
observes in the Monte Carlo history jumps between the Z2 states. Precisely both of these 
behaviours can be seen on the right hand side of Fig. [5l Until about bin number 400 the 
system was in the deconfined phase and then it changes into the confined phase. 

The location of the phase transition depends not only on but also on the ratio 
L/a/N^. In principle we would like to be in a situation where <^ L/ a m order to have 
a compact extra dimension. But the meaning of compactification will have to be qualified 
by looking at the results for the particle spectrum. 



4.2 Higgs and photon spectra 

In Fig. [6] and Fig. [7] the Higgs and photon masses for the ground state and the first excited 
state are shown as a function of (3 in units of 1/i? for the = 4, L/a = 8 and T/a = 96 
orbifold geometry. It was not possible to extract these masses for (3 < (3^ = 1.5975. In the 
confined phase the signal for the effective masses disappears immediately in noise. Right 
above the phase transition the signal for the particle masses is there. 

The first observation, looking at Fig. [6], concerns the Higgs ground state mass mh. For 
all P > Pc the Higgs mass is consistent with its value computed in 1-loop perturbation 
theory: for general gauge group G = SU{N) this is 

^ (4.1) 

where c = 3/{47t^)^/NC{3)C2{G)/2 and C2{G) = N {c = 0.1178 for SU{2)). Here we 
have taken the continuum result in [11] and replaced the dimensionful coupling by its 
lattice definition Eq. (12.31) . 



^ The expectation value of the basic Higgs field Eq. (|3.10p is actually constant for /3 < /?c 
The measurements for each simulations are blocked in 500 bins for the error analysis. 
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Figure 6: The ground and first excited state of tlie Higgs. Scan in (3 at fixed = 4. 
Masses are in units of 1/R. Tlie daslied line represents tlie l-loop result for the ground 
state. 
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Figure 7: The ground and first excited state of the photon. Scan in (3 at fixed = 4. 
Masses are in units of 1/R. 
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1.609 



1.80 



1.80 



1.65 



1.609 



1.65 



8 



8 



8 



12 



12 



12 



0.085(13) 
0.075(10) 
0.093(14) 
0.120(14) 
0.052(24) 
0.048(13) 



0.138(9) 

0.202(20) 

0.132(6) 

0.120(13) 

0.079(4) 

0.097(5) 



0.62(8) 

0.48(8) 

0.69(5) 

0.61(5) 

0.38(17) 

0.70(6) 



Table 1: Finite volume study of the spectrum at = 4. The excited state of the Higgs 
could not be determined at L/a = 12. 

The second observation, looking at Fig. [71 concerns the photon ground state mass 
m^. Contrary to the 1-loop prediction [31] the photon mass is non-zero for all f3 > f3c- 
The photon mass even increases as the phase transition is approached. This means that 
there is spontaneous symmetry breaking in the pure gauge theory. This is, to our best 
knowledge, the first non-perturbative evidence for the Higgs mechanism!"] originating from 
an extra dimension. 

The third observation, looking at both Fig. [6] and Fig. [TJ concerns the excited state 
masses for the Higgs mj^ and the photon m*. In perturbation theory, the first excited 
(Kaluza-Klein) states are expected to appear split from the ground states by (Am)i? = 1, 
the second with a mass splitting twice that, and so forth. In no range of f3 we see excited 
states at about 1 in units of 1/R. Close to the phase transition the excited states are 
separated from the ground state, especially in the case of the photon. Instead for larger 
(3 they get closer in mass to the ground states. This is an indication that at fixed N^, = 4 
the system behaves more like a compact system close to the phase transition rather than 
for large p. 

At this point one might worry about finite L effects in the particle masses. Especially 
so for the photon mass, since its non-zero value contradicts perturbation theory. In Tabled] 
we offer a comparison at several f3 values of the particle masses between L/a = 8 and 
L/a = 12. One can see that finite L effects are small and in most cases not significant. 
For sure there is no variation that could be explained with a behaviour m oc 1/L which 
is characteristic of finite volume effects. 

Spontaneous symmetry breaking in this context goes back to works by Hosotani [6,7]. 
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Figure 8: The static potential at (3 = 1.609 and = 4 between static charges on the 
boundary. The five-dimensional fits are clearly favoured. 

4.3 Static potential in the four-dimensional slices 

In this section we present simulation data for the static potential in the boundary slice, 
see Section 13.31 together with the results of various fits. Since we know from the results 
in Section 14.21 that the gauge boson associated with the U{1) gauge symmetry on the 
boundary is massive, the physically motivated fits are Yukawa potentials 

aV{r) = — Ci exp(— m^r)/(r/a) + Co in four dimensions , (4-2) 
aV{r) = —diKi{m^r) I ij- / a) + do in five dimensions , (4.3) 

where Cq, Ci and c?o, di are the fit parameters and the photon mass is the one mea- 
sured in the simulations. In Appendix [B] we derive the five-dimensional form of Yukawa 
potentials. It is interesting to compare the Yukawa fits with Coulomb fits 

aV{r) = —fi/{r/a)-\-fQ in four dimensions , (4.4) 
aV{r) = — ei/(r/a)^ + Co in five dimensions , (4-5) 

where /o, /i and cq, ei are the fits parameters. 

In Fig.[8]we present results from a simulation of the orbifold geometry T/a = 32, L/a = 
16 and = 4. The statistics is of 20'000 measurements of the potential. The photon 
mass is arrij = 0.108(7) (for the fits we neglect its error), measured in the simulations 
described in Section 14 .21 In the fits only the points r/a = 2, . . . , 8 are included. The 




data N^=4, n^=0 

5d Coulomb 
5d Yukawa 
4d Coulomb 
4d Yukawa 

7 8 



r/a 
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5d Coulomb 5d Yukawa 4d Coulomb 4d Yukawa 



4 0.4 0.6 14 10 

6 1.1 1.9 4.2 1.7 

Table 2: The values of x^/(degrees of freedom) for the various fits to the static potential 
at the boundary at /5 = 1.609. A comparison between = 4 and A^s = 6 is made. 
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Figure 9: The static potential at /3 = 1.609 and = 6 between static charges on 
the boundary. Besides the five- dimensional fits also a four-dimensional Yukawa fit can 
describe the data. 

per degree of freedom are listed in Table [2l The four- dimensional fits are excluded 
whereas the five-dimensional ones are very good. 

The data at = 4 might suffer from the fact that the ratio between the cut-off and 
the compactification scale is only R/a = N^/tt ~ 1.3. In Fig. [9] we present results from 
a simulation of the orbifold geometry T/a = 32, L/a = 16 and = 6. The statistics 
is of 9'000 measurements of the potential. The photon mass is am^ = 0.173(18), taken 
from a simulation at the same (3 and values but with T/a = 96 and L/a = 12. In 
Table [2] we can compare the per degree of freedom. We see that at = 6 together 
with the five-dimensional fits also the four-dimensional Yukawa form is a good fit. The 
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four- dimensional Coulomb fit is instead excluded. It should be noted that the errors on 
the potential data are comparable at both values of Nr,. 

We also measured the static potential in the middle slice at = 2. The results of 
the fits confirm the conclusions from the fits in the boundary slice at = 4, that is 
the four-dimensional forms are disfavoured. The values are now about 3 for the four- 
dimensional forms and 0.3 for the five- dimensional forms. At A^s = 6 the statistical errors 
are too large and all the fits are equally good. As was discussed in Section 13.31 we have 
to improve the measurements of the static potential, especially for the middle slice. It 
is important to have precise data for the latter, since any difference with respect to the 
potential on the boundary might be a hint of localization effects. 



4.4 Discussion 



A certainly surprising feature of these models emerges when we look at the mass spectrum 
itself, measured in i?-units. We have seen in Fig. [6] that when the compactification scale 
is lower than the cut-off scale by a fixed gap (i.e. for fixed), the Higgs ground state 
mass is consistent with its 1-loop perturbative value not only at large values of P but also 
in the definitely non-perturbative region when /3 is close to its critical value at the phase 
transition. Moreover also the photon mass has a mild dependence on /3, see Fig. [71 Since 
we do not have at present any powerful analytical tools to explain these observations, we 
will keep our discussion at a qualitative level. 

We would like to argue that one way to interpret this peculiar phenomenon is by a very 
mild cut-off dependence of the five-dimensional bare coupling in the non-perturbative 
regimi^. The generic situation observed in our simulations, as far as the spectrum is 
concerned, can be translated in an effective field theory language by saying that the 
theory possesses a region in its parameter space such that any operator of dimension 5+p 
appearing in the effective action and effective operator^ scales as 

with Z{gQ,g4) a very slowly varying function of the dimensionless bare couplings go = 
g^VX and g4 = g^/V tcR^ at least as long as the ratio A/(l/R) is kept fixed. In the orbifold 
theory, taking the Higgs mass as an example, we recall that a boundary counterterm is 
non-perturbatively excluded and thus all possible corrections come from either pure bulk 



^^Note that this is not in contradiction with the perturbative, 1-loop expectation which wants the 
five-dimensional bare couphng to have a mild cut-ofF dependence for energy scales E ~ \/R where the 
theory is truly five-dimensional (see eq. (3.13) of [53]). 

^■^ This is Symanzik's analysis of cut-ofF effects. 
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effects or bulk effects with boundary insertions, both of which descend directly from the 
circle. 

In an attempt to predict the explicit cut-off dependence in Eq. (14. 6 p , we note that naive 
dimensional analysis tells us that as j3 decreases with fixed, the cut-off increases, see 
Eq. (12. 3p . The compactification scale is l/i? and a wide separation from the cut-off scale 
requires A >> 1/R. Increasing (3 while keeping the gap between the compactification 
and cut-off scales fixed, would require decreasing A (at fixed g^) and therefore an increase 
of which drives the fifth dimension to its decompactification limit. A general lesson 
then is that for fixed A^s, moving towards the large (3 regime is expected to enhance the 
cut-off effects (appearing as -E/ A at low energies E in the sense of an effective action) and 
decompactify the theory, whereas moving in the opposite direction, i.e. towards smaller 
/?, is expected to suppress the cut-off effects and drive the system into a compactified but 
non-perturbative regime. Eventually the phase transition is reached at the critical value 
of /3 = (3ci where the cut-off reaches its maximal value. There is a possible caveat though 
in this argument: we have implicitly assumed that it makes sense to vary the cut-off while 
keeping the dimensionful coupling g^ fixed for all values Pc < P < oo. 

If we look at the spectrum of the excited states of the Higgs and the photon in Fig. [6] 
and Fig. U\ we see that as (3 grows their masses decrease, thus deviating more and more 
from their expected perturbative value as Kaluza-Klein states of mass ~ This fact 
could be explained with our argument that cut-off effects increase at large (3. We cannot 
at the moment explain why these cut-off effects would affect the excited states but not 
the ground states. 

Next, let us see what happens when we start changing while keeping j3 fixed. The 
cut-off A depends only on f3 and is also fixed. Using Eq. (12.21) we obtain that 

6Nr, SR 



N. R 



(4.7) 



which implies that increasing Nr, amounts to increasing also R. Thus, in a compactified 
scenario we expect to be able to approach an effective dimensional reduction by decreasing 
Nr,. It is interesting to note that according to our potential data we observe an opposite 
effect: the fit to a four dimensional Yukawa law becomes much better when we increase 
A^5 from 4 to 6 at fixed (3, namely when we decompactify the extra dimension. Therefore, 
to the extent that this result can be considered as a firm physical property of the system 
in this region of the parameter space, we are lead to the conclusion that the fact that 
we start observing an effective dimensional reduction at A^s = 6 is more likely to be a 
consequence of a localization mechanism rather than an effect of compactification. It is 
possible that the region that corresponds to dimensional reduction from compactification 
is located at much smaller which would however require an anisotropic lattice to be 
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probed. 

Finally, regarding the applicability of the vicinity of the phase transition in model 
building, we would like to point out that even though it clearly corresponds to a strong 
coupling regime from the point of view of five dimensions, viewed from the point of view of 
the four-dimensional effective theory with an effective coupling defined as in Eq. (12. 4p . it 
could correspond to a weakly coupled regime, if for example large renormalization factors 
change Eq. fl2^ . 

5 Conclusions 

The first results for the spectrum of the orbifold theory simulated on a lattice with an 
extra dimension of the size of = 4 lattice spacings, show that there is spontaneous 
symmetry breaking, which manifests in a massive gauge boson associated with the U{1) 
boundary gauge group. This result was unexpected from computations at 1-loop in per- 
turbation theory, where the photon remains massless. Moreover, the Higgs mass is also 
measured in a large range of gauge coupling and it is always close to its perturbative 
value. The excited state masses for the Higgs and the photon are not at the scale of the 
inverse compactification radius, where they are expected to be in a scenario of dimensional 
reduction like in finite temperature field theory. But close to the orbifold phase transition 
the mass splitting between the ground states and the first excited states increases. Data 
for the static potential strongly suggest five- dimensional potential forms, both when the 
static potential is measured with the charges on the boundary slice or in the middle slice 
along the extra dimension. We have also presented potential data at A^s = 6. They 
indicate that a four- dimensional Yukawa form, with the vector boson mass equal to the 
measured photon mass, is a good fit to the potential on the boundary together with the 
five-dimensional forms. 

The results we have so far give a consistent picture and they motivate for further work 
to explore the phase diagram of the orbifold theory. Despite its non-renormalizability, 
the theory makes finite predictions for the mass spectrum and we would like to under- 
stand better how to change the lattice parameters to study the scaling properties. This 
also requires technical improvements, for the basis of Higgs operators and for the static 
potential. 
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A Higgs potential from extra dimension(s) 

In this Appendix we review the calculation of 1-loop potentials in orbifold gauge theories. 
This material is already well known in the hep-ph community but since it is perhaps less 
familiar in the lattice community, for completeness we reproduce here in detail the pertur- 
bative arguments for the existence or not of spontaneous symmetry breaking, reproducing 
essentially some of the the results of [31]. 

Consider a (massive) free field theory for a one-component real scalar field in D- 
dimensions with action S. After a Euclidean rotation, the path intergral that defines the 
vacuum energy F is [8]: 

e-r = f [Dcf>] e-^^ ^ , □ = . (A.l) 

The mass (M) dependence of the above can be extracted by using the identity 



lo, 

We obtain 

r = -log 



-tr (e-*-^) . (A.2) 



1 



2(2ir)D7„ t 



D+2 
2 



Vdet [□ + M2] 
For a fermionic free field we have 

e-r^detp), D^T^a^ + M. (A.4) 

The nonzero eigenvalues of the Dirac operator come in complex conjugate pairs ±ia + M, 
hence 

det(D) = ydet(DDt) = ^detp + M^] . (A.5) 
We can therefore summarize the contributions to the effective action, setting D = 4, as 

where Fj is for bosonic and 1 for fermionic degrees of freedom of mass M/. 
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Let us apply the above for the case where there are d extra compact toroidal direc- 
tions. For the four- dimensional theory, after dimensional reduction, i.e. Kaluza-Klein 
decomposition, we obtain the mass formula 

M,^ = m? + X:(^)', (A.7) 

where 



• Ri are the radii of the d circles of the torus. 

• {rii}, i = 1, . . . ,d are integers. 

• mf is a 4 -|- d-dimensional mass which remains in the Ri — > oo limit. This is zero 
for a pure gauge theory. 

• The shifts af originate from the possible failure for periodicity of the i + d dimen- 
sional field $/: 

^i{x'',y' + 2nhRi) ^ e^'^^''^''i^i{x'',y'), (A.8) 
where yi are the coordinates of the circles. 

According to the above, for the torus we have the expression for the potential 

ySf = -ET,(-f^J d/ie-^'^. (A.9) 

By commuting the integral with the sum over rii and then performing a Poisson resum- 
mation using the formula 

g-TT nT'An+2i7rb'rn ^ ^ ^-tt (m-b)'rA-l (m-b) 

{ui} {mi} 

where A is a d x d (invertible) matrix and n and m are d dimensional KK number vectors, 
we find 

I 327r 2 ^^^^ Jo 

The terms with nj = give rise to a divergent contribution to the vacum energy. They 
represent a contribution to the cosmological constant and can be neglected for the present 
discussion. For all other non-vanishing vectors n we perform the integral explicitly. This 
leads to the finite result 

I 327r 2 \^ , y n^o(Ei^f^f) ' 
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A.l 5D SU{2) gauge theory on S^/Z2 



For a 5D theory compactified on the 1 dimensional torus, the circle, we have d = 1. Then 



r(i) 



2Trina 



T/51 ^ _ \^ ( _\Fi p \^ J: ^ '3 "sp 



n=l 



COS {2Tma^ 



(A. 13) 



For a pure gauge theory we have in addition that Fj = 0. We then finally obtain for the 
vacuum energy 



00 



_ ^^ C0S(2W; 

A n=l 

where the index / has been changed to A, indicating the adjoint representation of some 
gauge group. This potential may have a minimum that leads to an extra dimensional 
version of the Higgs mechanism. 

There are different cases where a failure of periodicity can happen and a shift in the 
KK formula is generated. One particularly interesting example is the presence of a Wilson 
line in an extra dimensional gauge theory 



a'' = q''9, / TT^A^, (A. 15) 

where is the internal component of the gauge field with gauge coupling and qa the 
charge of the Ath field under the corresponding generator. To be more concrete, let us 
look at an example. Consider a 5D SU (2) gauge theory compactified on the circle. In 
order to compute the vacuum energy, according to eq. ( 1A.14I) the only quantity we need 
to compute are the a^. To compute the we have to compute the eigenvalues of the 
mass-squared operator [31] 



— DmDm = —D^D^ — D^D^ (in Euclidean space) (A. 16) 

in the vacuum, with Dm acting on fields F = F^T^ in the adjoint representation as 

DmF = OmF+IBm^F], Bm = {Am). (A.17) 

Bm is a constant background field given by the vacuum expectation value (vev) of Am- 

After orbifolding, which breaks SU{2) — > U{1) the only fields that can take a vev 
are Al{x^) and Aflx^). The unbroken U{1) invariance can be used to rotate the vev's in 
such a way that we have 

Bl = H, Bl = Bl = 0, B, = -ig^B^T^ . (A. 18) 
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The non-trivial part comes from the D5D5 part of the operator acting on dependent 
parts of the gauge fields. Writing out this explicitly in the vacuum acting on an adjoint 
field F, 

D5D5F = d5d5F + 2[B,,d,F] + [B,,[B5,F]], (A.19) 

gives the operator 

/d,d5 \ 

{D,D,)^j, = ^5^5 - glHH -2g,Hd, . (A.20) 
V 2g,Hd, d,d,-glHH/ 

Since this operator does not mix different KK modes we can diagonalize it separately for 
each level n. 

• Eigenvalues of A^. 

The matrix elements of the —D^D^ operator 

/•27rR 

{f\-D,D,\g) = / dx" r{x'){-D,D,)g{x') (A.21) 
can be easily obtained by using the basis of orthonormal functions 

1 Tl 

cos for A,B^?> (A.22) 



smce is even under the orbifold, and 



1 n 

sin for A, B = 1, 2 (A.23) 



since ^4^'^ are odd under the orbifold. The matrix we obtain for n 7^ is 

v? {\n\+af {\n\-af ^4) 


















B? ii2 


9 an 




u 


r) an 

^i?2 


R? il2 



i?2 ' R2 

where 

a^g^HR (A.25) 

and on the right we have given the eigenvalues. For n = 0, the ^4^'^ do not have 
zero modes but A^^ does have a zero mode whose eigenvalue is 

^. (A.26) 
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• Eigenvalues of ghosts. 

We notice that the ghosts have the same Z2 parity assignement as the gauge fields 
Ay, therefore the effect of the ghosts is to just reduce the degree of polarization as 
4-^2. 

• Eigenvalues of A^. 

Since the parities of the A^ are the opposite from those of A^, we will obtain 

g (A.27) 

and 

.A(n^o) (1^1 +«)^ (1^1 -Q^)^ ofi^ 

By comparing with Eq. (1A.7I) we see that non-zero shifts are given by ±a in Eq. ( 1A.24I) 
and Eq. (lA.28p . As we have argued, only the non-zero modes are relevant for the vacuum 
energy. We count 2 + 1 (2 from physical degrees of polarization of A^, 1 from A^ 
eigenvalues in + a)^ jR? (with a} = a) and an equal number of eigenvalues {n — a^/R"^ 
(with = —a). Furthermore, as can be seen from eq. flA.14p the = —a contribution 
is the same as the = +a contribution. We then obtain for the vacum energy the final 
result 



vSASU{2)) = -^^4:4 y h^^^A^ . (A.29) 

n=l 

This potential could, in principle break the remaining U{1) down to nothing on the branes. 
Plotting the potential one can see that it has two degenerate minima at a = and a = 1, 
none of which breaks U{1). 

A similar computation gives for the SU{3) — > SU{2) x f/(l) model the potential 

KffiSUiS)) = E (2™) + 2 cos {nna)] . (A.30) 

71=1 

The latter can be seen by noticing that the eigenvalues for the adjoint of SU (3) are 

(n±a? (n ± f )2 ^ 

2X-, ^ 2x^. (A.31) 

Again, the minimum is at a = which does not break SU{2) x f/(l). 

In fact, in some cases a useful criterion to see if it is possible to break the symmetry by 
the Hosotani mechanism is to use the statement proven in [11]: The Hosotani mechanism 
does not reduce the rank of Ti. if the symmetry breaking global minimum is at a = 1. To 
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show this statement, we can compute the Wilson hne due to the vev H — ^ of a scalar 
along the direction: 

It is straightforward to show that exp(— 27riT^) is a diagonal matrix. In the case of SU (2): 
exp(— zTrcTi) = — 12 = exp(— i7r(T2). For SU{N), the non-diagonal generators are obtained 
embedding ai^2- Thus, we always have that 

[{W),H,] = Q, (A.33) 

where Hi are the generators corresponding to the Cartan subalgebra of G, i-e. that the 
Wilson loop commutes with at least those generators and therefore it leaves at least a 
C/(l)i X • • • X U{l)rank{H) unbrokcn. Prom this it is clear that in the SU{2) model, C/(l) 
can not break further with a — 1. For the SU{3) model a — 1 could break SU{2) x U{1) 
down to U{1) X U{1) by the Hosotani mechanism. In both cases, it is not clear what 
a ^ 1 would do. 



B The Yukawa potential in 5D 

The potential is given by the integral 

/j4 iq-x 



where \x\ — r, q — \q\q', \q'\ — 1 and 



q'l — cos 01 

^2 = sin (pi cos 02 

^3 = sin 01 sin 02 cos 9 

^4 = sin 01 sin 02 sin 9 (B.2) 



where < 0^ < tt, < ^ < 2n. Then the potential can be written as 

-j^ POO p2tt ptt pit ^iq-x 

^i^)^77^ / d\q\\qf d9 / c?0i / c?02 sin^i sin 02^— — (B.3) 
l^'^j Jo Jo Jo Jo q +m 

Choosing coordinates such that x = r ■ ci with ei a unit vector, we have that q ■ x — 
\q\r cos 01 and then 





111 


3 




2 _|_ 



/ de #2 sin 02 / d0isin2 016^1*1^ ""^"^i. (B.4) 
Jo Jo Jo 
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Performing the angular integrals we obtain 

Finally, changing variables as y = \q\r one has the left over radial integral 

^(^) = r dy^-^^My): (b.6) 

which can be done, yielding the result 

- ^^^(mr)Xi(mr). (B.7) 
As r — > 0, the Bessel function {mr)Ki{rnr) — > 1 so we have at short distances 

- ^ : V{x) = (B.8) 
For large r on the other hand we have 
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